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I complete in this paper the solution of the Bukhvostov Lipatov model by computing 
the physical excitations and their factorized S matrix. I also explain the origin of the 
paradoxes which led in recent years to the suspicion that the model may not be integrable. 
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In a famous paper of 1981, Bukhvostov and Lipatov (BL) |IJ presented pieces of the 
solution to a very interesting quantum field theory model, whose bare Lagrangian reads 

C = -01 (i$ - m) -01 + ^2 (i$ -m)ifj 2 + 9b$\l '^Y^il^i- (1) 

Except for the massive Thirring model, most known integrable models with several fermion 
species do not contain an explicit mass term: the Gross Neveu models, SU{2) and U(l) 
Thirring model etc rather all exhibit spontaneous mass generation. The model ([]]) is quite 
different in nature and it has not been clear up to that day what sort of "family" of 
integrable theories it belongs to. 

In their original paper, BL succeeded in diagonalizing the model (|1]) with the coordi- 
nate Bethe ansatz method. They did work out only half of the solution however, building 
the ground state, but stopping short of discussing the physical excitations. Although there 
are well established methods to do so in principle 0, 0, in practice, the bethe equations 
written by BL are complex, leading to a bewildering array of excitations, where it is hard 
to identify fundemental particles, bound states, and "pseudo-particles" arising from non 
diagonal scattering. Certainly, the lack of a quantum group understanding of the equations 
in JI] did not help in that endeavour. 

In the last few years, a growing suspicion has mounted that maybe the BL model 
was not solvable after all. When naively bosonizing ([]]), one obtains a double sine-Gordon 
model which indeed appears to be not integrable classically ||. Worse, this bosonized 
model also appears to be not integrable (this aspect has again been emphasized very 
recently in ||) within the framework || of conformal perturbation theory! My purpose in 
this note is to show the simple way out of these difficulties, and to finish up the solution 
of the BL model by computing the physical excitations and their factorized S matrix. 

To start, I recall the Bethe ansatz solution of [0J. Writing the energy 

E = m sinh yj , 

j 

the equations quantizing allowed momenta are 

/ • t ■ i \ TT sinh fe -*r + ij) 
exp {-imL smhyj) =11 ——, 

-A- sinh (z s - yj + ig) J-r sinh (z s - z r + ig) 
sinh (z s - yj - ig) ^1 s i n h _ Zr - ig) 

1 



(2) 



(3) 



In the latter equations, I have sligtly changed notations compared with [fl]. I have intro- 
duced the coupling g, whose exact relation to the bare coupling in ([D depends on the 
regularization used in solving the coordinate Bethe ansatz equations (it is, quite confus- 
ingly, denoted by the same symbol in |1| however!). Compared with equations (81,82) in 
I have also switched the sign of the coupling constant, and will moreover restrict to the 
case g > here (that is, g < in notations of [0). I have set their v r = z r — y. Finally, I 
have chosen antiperiodic boundary conditions for the fermions. If I call Ni the number of 
fermions of type i, then N2 = I and iVj = n — I in (|3|). 

While BL used a sharp cut-off regularization, I would like to proceed slightly differ- 
ently, and introduce the cut-off A at this early stage, as was done for the Thirring model 
in J7j. One has sometimes to be careful with the choice of cut-off: for instance, it leads to 
drastic differences in the physical properties of the Thirring model in the repulsive regime 
0,0. I have however checked in this particular case that whatever procedure gives iden- 
tical results; the reason why I want to use the smooth cut-off of is to facilitate later 
comparison with lattice models. I will thus replace the left hand side of the first equation 
in (H) as 



exp (—imL sinh y) 



I L 1 2 d 

sinh 2 (yj — A + z|) sinh ^ yyj + A + z|) 
sinh § ( Vj - A - q ) sinh \ (y 3 - A - %\ ) 



(4) 



When A is large, a (having the units of length) is small and y « A, this reproduces the 
previous term, together with the correspondence 

m = 2 sm 777 (5) 

while divergences are smoothly cut-off at large value of y. I similarly will use for the energy 
the derivative of the momentum read off from (f|), that is 



v, 1 • 9 

m cosh Vj — > - sin - 

yj 2 2 



1 1 

+ 



cosh(?/j + A) — cos I cosh(yj — A) — cos | 



(6) 



I can now proceed and study the physics encoded in these equations. As pointed 
out in |l[], it is easy to check that the ground state is obtained by filling up a sea of yj 
antistrings, together with a sea of real (1-string) z r . A first possible strategy from then on 
is to study the possible physical excitations obtained by making holes, adding other types 
of strings etc, and try to extract their S matrix. This turns out to be a rather confusing 
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task however, for a reason that we will understand easily later: whatever g, the scattering 
is not diagonal, and there is a complex spectrum both of bound states and pseudo particles. 

To make progress, I rather will use the approach pioneered by Kliimper and Pearce , 
and, independently, by Destri and de Vega, jlQ|] , and widely used since to tackle theories 
with complicated scattering |TI |. After the usual manipulations, the KPDDV equations 
read 

f{y) =iL2M cos ^ sinh ?/ + 2z J dy'^n{y - y')Im In (l + e" /(?/ - l0) ) 
+ 2i [ dz$ 12 (y-z)lm In (l + e"^-* )) 

r , (7) 

g{z) =iLM sinh z + 2i / dz'^niz - z')lm In (l + e ~ 9{z '~ i0) 

+ 2i J dy$i 2 (z - y)Im In (l + e ~ f(v ~ l0 ^ 
where I have set M = — ^ir^. The energy then reads 



L 

E = — 

7T 



2Mcosy J cfysinhylm In (l + e _/ (y_i0) ) + M J dzsmkz Im In (l + e ~ 9(z ~ i0) 

(8) 

In (|7|) and (|8|), the integrals are running from — oo to oo. In these equations, the kernels 



are given by, setting g = , t a real number, 

sinh ^ ^ x 
$n = — 

1 1 - (t+2)x 



sinh 

$ 22 = J^f x ( 9 ) 

sinh ^ ^ x sinh ( t+2 ^ x 



2 2 

sinh y 



(t+2)x 



sinh , 



where we have introduced the Fourier transform f(x) = ^- J f(y)e' ltxy ^ 7r . 

Two things can be rigorously deduced from these equations. The first is that the UV 
limit (to — > 0) of the theory has central charge c = 2, as expected from (P. The second 
is that the physical mass is simply proportional to the bare mass, not a power of it: this 
means that in the physical (renormalized) theory the BL equations are describing (I will 
get back to this issue later) , the operator perturbing the UV fixed point must have scaling 
dimensions x = 1, irrespective of g. 

Besides, extracting the scattering theory from the KPDDV equations is still a matter 
of guess work. So far, these equations have had the very simple feature that they contain 
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only the most "basic" ingredients of the theory: the fundamental particles and their S 
matrix. In the sine-Gordon case for instance the right hand side would involve only 
one type of terms (only one distribution function) , with kernel $ = i In S++ , S++ the 
soliton-soliton scattering matrix, y the rapidity. 

In the present case, I claim we can interpret the equations as follows. First, I observe 
that, writing 

sinh 2 x sinh x sinh x 

sinh (*=2)£ sinh ~ 2 cosh % sinh ~ 2 cosh 5 sinh ' 



2 2 ^ 2 2 



we can identify the $22 terms in the KPDDV equations as $22 = |r In 
I have introduced the two parameters (recall g — ^f-) 

t-2 



Here, 



(3( = 4tt- 



t- 1 



52 - ^ + 2 



(10) 



t + r 

while by S@ I denote the S matrix for an ordinary sine- Gordon theory whose parameter is 

4tt ' 



(3 [|12| (that is, the dimension of the perturbing operator is x = §-) 



S+ + (y) = exp 



<i£ . 2y£/i, sinh((U — 1)£ 
2£ tv sinh £ cosh 



8tt-/3 2 

^=^2— (11) 



This leads to my basic guess: the physical theory is made up of four particles (kinks) of 
mass M carrying a pair of quantum numbers Qi, Q2 = ±1 (how these are related to the 
original problem will be discussed soon), and which scatter with the S matrix 

S = S&®S&. (12) 

Notice that p\ < 4tc while /3| > 4tt for t e [2, 00). Therefore, while the second S matrix in 
C pL2f) is in the repulsive regime, the first one is in the attractive regime, and therefore will 
exhibit bound states. In the KPDDV equations, neutral bound states do not show up. Here 
however, because there are two charges, it is reasonable that the "basic charged" bound 
states should also appear. In fact, one easily checks that the second mass in our equations, 
2Mcos ^2, is precisely the mass for the first bound state in a SG theory with scattering 
S^ 1 . Moroever, one can also check that the kernels $n and $12 are the exact kernels one 
would obtain when scattering one of our bound states with either another bound state, or 
a basic kink. In doing this check, one should not forget that, although S^ 2 has no bound 
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state, it of course does contribute to the overall scattering of bound states. For instance, 
$12 arises from the scattering matrix Sf^_(y)S+ + (y — ^) (y + where Sb+ is the 
soliton one- breather S matrix in the usual sine-Gordon model with parameter /3i. 

The claim (|12|) therefore appears at least reasonable from that perspective. What I 
have done next is go backwards, and checked it carefully against the Bethe equations by 
using the more traditional method of identifying the basic excitations and computing their 
scattering. Equation (|T2]) turns out to be perfectly confirmed. I found out in particular 
that making a hole in the z distribution produces a fundamental kink, while making a hole 
in the y distribution produces a fundamental bound state (observe this is true provided 
t < oo, that is g > 0. The free limit is singular from the point of view of the Bethe 
equations, which does not help in analyzing them). Other bound states are obtained by 
complex distributions of roots (in particular "strings over strings" ) which do not seem too 
interesting to discuss here. An additional result (which could also have been obtained by 
adding magnetic fields in the KPDDV equations) is the charge of the fundamental kinks 
in terms of the original charges: one has 

Qi = ^(iVi + iV 2 ) 

2 J (13) 
Q2 = % L {N 1 -N 2 ). 

Z7T 

Here, I have introduced still other parameters, for a reason that will hopefully become 
clear soon: 

Pi _ hi _t-2 

2tt 8tt - ft t 

ft ^ hi = t + 2 (14) 
2?r 8n-hi 1 

Notice that as t — > co, Qi ^A^i+ A^2, Q2 —> -^1 — N2, so the basic particles coincide with 
the four elementary fermions in that limit. 

At this stage, the alert reader will have no doubt recognized that the scattering theory 
we have extracted from the BL equations is nothing but the scattering theory for the double 
sine-Gordon model [0 (a particular case of a general model studied by Fateev |i~4|), whose 
renormalized Lagrangian reads, after bosonizaton 

C = -d^fa + \d^4>2 + Acos cos/? 2 </> 2 , (15) 



with 

/3 1 2 +/3 2 2 =4tt. (16) 

The notations are of course chosen so that (|16D matches (|I4]) . In particular, the conditions 
( |16|) means that the dimension of the perturbing operator in ([15|) is x = 1. 

In fact, this result is not so surprising. It is easy to check, at least if one requires the 
existence of a conserved quantity of spin 3, that the only non trivial manifolds where the 
double sine-Gordon model is quantum integrable are given by and, maybe, the other 
manifold f3f + (3% = 4n. 

From that perspective however, the naively bosonized theory associated with ([]]) cor- 
responds to 4j + 4? = -, and looks completly baffling! What happens is, I think, quite 
simple. In the coordinate Bethe ansatz, one always deals with a bare Lagrangian, that is 
then regularized by using a particular prescription: making sense of terms like 5(x)sign(x) 
when one solves the Bethe equations, and introducing a cut-off in the rapidity integrals. 
There is of course no reason why the resulting large distances properties should be de- 
scribed by a renormalized theory whose parameters are the same as the bare ones! This 
well known fact is hammered home by the example of the ordinary Thirring model (with 
four fermion coupling gj<): Korepin for instance [|J found that (3 2 = A{rc — gj>), (Bergknoff 
and Thacker ]T5[ still have a different result) which differs from Coleman's [|IB| famous 
correspondence (3 2 = ^ 9T . None of this is surprising in the least of course, and it did 
not matter very much so far, because the models one was dealing with were always in- 
tegrable anyway; it was just a matter of knowing which particular point in one language 
corresponds to which particular point in the other. 

Things are very different here, since the models of interest are integrable only in a 
subset of the whole parameter space: Lagrangians have to be specified much more carefully, 
and the integrable theory might look quite different depending on which point of view one 
is adopting. In bosonization as well as in conformal perturbation theory, one usually deals 
with renormalized theories. There is thus no reason why, by naively bosonizing the bare 
Lagrangian of Bukhvostov Lipatov and interpreting it at face value for a renormalized 
Lagrangian, one should find a theory that is integrable in conformal perturbation theory. 
In other words, the bare Lagrangian they wrote, together with the regularization they used, 
define an integrable quantum field theory, and in that sense, the Bukhvostov Lipatov model 
is integrable 0. But it was misleading of the authors in to proceed with bosonization, 



1 Classical integr ability, as checked in S, also follows. 
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and the final double sine-Gordon model they wrote down, with -hz + 4k = —, cannot be 
expected to be integrable. 

The only proper way to proceed, once faced with (]!]), is to identify the scattering 
theory by studying the bare and physical Bethe ansatz equations. Once the S matrix 
( |12D is obtained, one can for instance observe that it has affine quantum group symmetry 
sl qi (2) <S> sl q2 (2) [HJ; this, together with the fact that the dimension of the perturbing 



operator is x = 1, leads unambiguously to ( |15|) with (^6|). After refermionization, (jl 



reads as (Q) with the additional appearance of terms (^7^^) . These terms come with 
a coefficient of order g% at small gr& - as in the Thirring model, renormalization effects are 
seen only at higher orders, and at leading order in the bare coupling constant all models 
are equivalent. 

To make things more concrete and somewhat more rigorous, I would like to finally 
point out that the problem I have been discussing can be studied quite explicitely with 
a lattice model regularization. As discovered in [|17|j , the following hamiltonian, obtained 
from a twisted Osp q {2/2)^ R matrix [fT8| , is exactly solvable: 

H = ^ ( C l> C i+ 1 > £T + CC ) i 1 - n 3-° ~ n 3+l -a ~ aV l ( n 3-<y - n j + l,-a)) 

+ V 2 J2 (c},+ct _c i+ i,_c i+ i >+ - ct )+ c] +1 _c j+li+Cj> - + ccj 

3 

+ V 2 ^2 i n 3,+ n 0- + n 3+l,+ n 3+l,- + n 3,+ n j+h- + n j,- n 3+l,+ ~ n 3 ~ n 3 + l + X ) > 

(17) 

where n JiCr = cj a Cj >(T , V\ = sin 7, ]/% = cos 7, q = e 11 . By using the well known techniques 
(see eg JHJ) to take the continuum limit of this model at small Vi, V2 (both negative) and 
half filling, and keeping only the relevant or marginal terms, I have found that ([H]) exactly 
gives rise to ([j]), with gj, cx — V2 and m = (the complicated fine tuning in (|PT|) cancels 
out the terms that would induce a gap in the similar looking Hubbard model). Of course, 
there is an infinity of additional irrelevant couplings, that will give rise to renormalized 
coupling constants: this is very similar to what happens in the XXZ model say, but here, 
this renormalization changes the form of the naive interaction quite a bit. By using the 
Bethe ansatz equations written in HPT] , I have checked that the CFT associated with (|T7| ) 
corresponds to the A —> limit of ([HI). It is in fact possible to put a mass term in the 
lattice model too by using an inhomogeneous distribution of spectral parameters as in [2D ]; 



the bare hamiltonian looks then as (]!]) , while the Bethe equations are identical with those 



I used before (A being then, as in |2C], a measure of the inhomogeneity, and a the lattice 
spacing) . 

In conclusion, it is a bit disappointing to realize that we have only one integrable 
manifold in the double sine-Gordon model @, the appealing but mysterious one hinted at 
in [|l| finally coinciding, after proper analysis, with the one in ||13| , ||14| . On the other hand, 



I hope that this discussion will lead to further progress in understanding theories with 
several bosons. As far as the 0(3) sigma model is concerned, it seems that the Bukhvostov 
Lipatov approach was almost right after all; according to recent work of Al. Zamolodchikov 
PHI , the proper theory describing the instantons anti-instantons interaction differs from 
(|15D , Qilf) , by the simple replacement /?2 — > ifii- 
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2 The are growing indications that the manifold /3? + /3f = & 71 is n °t integrable; note that the 
conformal perturbation theory argument is rather weak in that case, due to the operator having 
dimension one. 
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